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Abstract: We study the scale dependence of the non-linearity parameters /nl and 
#nl in curvaton models with self-interactions. We show that the spectral indices 
rif NL = dln|/ NL |/dln k and n gNL = dln|g NL |/dln/c can take values much greater than 
the slow-roll parameters and the spectral index of the power spectrum. This means 
that the scale-dependence of the bi and trispectrum could be easily observable in 
this scenario with Planck, which would lead to tight additional constraints on the 
model. Inspite of the highly non-trivial behaviour of /nl and ^nl m the curvaton 
models with self- interactions, we find that the model can be falsified if (7nl(^) is also 
observed. 
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1. Introduction 

Non-Gaussianity of the primordial perturbations can efficiently discriminate between 
different models of inflation. It is by now well known that both the strength and sta- 
tistical properties of primordial non-Gaussianities depend crucially on the details of 
the inflationary model. While conventional slow roll models of inflation with canon- 
ical dynamics typically predict negligible non-Gaussianity, non- minimal construc- 
tions may generate observable non-Gaussianity. For a selection of recent reviews see 
[0, 0, |3], f|, |5[ . Examples of such scenarios include models with non-canonical kinetic 
terms or a breakdown of slow roll dynamics and models where the primordial per- 
turbations are generated at the end of inflation, like in modulated reheating || [7]], 
or after the end of inflation, like in the curvaton scenario [|8]. |9| . 

The simplest type of non-Gaussianity is the so called local form. The local 
Ansatz for the primordial curvature perturbation reads 

Ck = C k G + ^/nl(C g * C G ) k + ^nl(C g * C G * C G ) k , (i.i) 

where ( G is a Gaussian field and the non-linearity parameters /nl and #nl are con- 
stants. The star denotes a convolution. In [ITUl O it was shown that non-linearities 
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of the field equations in general give rise to a mild scale dependence of /nl and c/nl 
even in models with canonical slow-roll dynamics during inflation, which typically 
have been analyzed using the local Ansatz. Deviations from the local Ansatz were 
analyzed previously from an observational point of view in [12, H| and more recently 



in |14|, |T3J. The scale dependence of /nl and #nl can be described by the quantities 

_ dln|/ NL | _ dln|g NL | 

n/NL " dlnfc ' n9 ^~ dlnfc ' [ ' 

and in models with slow roll dynamics during inflation the typical magnitude of 
rif NL and n gNL is set by slow roll parameters [11]. Hence the local Ansatz (|1.1|) gets 



replaced by a quasi-local form with /nl and #nl being weakly fc-dependent functions. 
Although the scale dependence is typically weak (for an exception see [0), it could 
be an observable effect which makes this topic very interesting. This topic is also 
interesting in models of non-local Gaussianity, see e.g. [|H], [L8|, [19], p0| . 

In this work we compute the scale dependency of the non-linearity parameters 
in self-interacting curvaton models. The scale dependence in the limit of weak in- 
teractions, where the self-interaction does not dominate over the quadratic part in 
the curvaton potential, was considered already in ||21|| . Here we analyze also the 



self-interaction dominated regime, which turns out to be particularly interesting. 
In this regime the non-linearity parameters /nl and #nl depend sensitively on the 
curvaton value at the time of inflation |22, ^3], |24|, ^5], p6fl . In particular, they os- 



cillate around the naive estimates |/nl| ~ r dec an d |<?nl| ~ r^ c , where r^ec denotes 
the curvaton contribution to the total energy density at the decay time. We find 
that the non-linearity parameters become strongly scale-dependent in the regions 
| /nl I ^dec, Ipnl I "C r^ c . For small values of r^ec the non-Gaussian effects in these 
regions are at the observable level and the scale dependence will be a detectable fea- 
ture of this model in the near future. We also comment on how the results could be 
generalized to any models where the primordial perturbation arises from a component 
which is sub dominant during inflation. 

The paper is organized as follows: in Section |2] we review the self-interacting 
curvaton scenario to the extent needed in our analysis. In Section [5] we derive expres- 
sions for the scale dependent non- Gaussianity in the curvaton scenario and discuss 
their generic features. In Section | we discuss the case of quartic self-interactions, 
deriving both analytical and numerical results. In Section ^ we discuss the scale 
dependence in curvaton models with non-renormalizable self-interactions, while in 



Section |57y we study the regime of very large self interactions and find an analytical 
approximation. Finally, we present our conclusions in Section ||. We use the units 
M P = (8ttG)~ 1 ' 2 = 1 throughout the paper. 

2. Self-interacting curvaton model 

We consider self-interacting curvaton models where the curvaton potential is given 
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by 



V = \m 2 a 2 + Xa n , (2.1) 



where n = 4,6, 8. 

While the curvaton should oscillate before decaying, and the quadratic part will 
typically be dominant at this stage due to the small field value, the self-interacting 
part may play an important role at earlier stages, crucially affecting the predictions 



of the curvaton scenario p2l E3L E4L EE 26, B7L EE E9, p(J, 0' 32]. We assume 



the primordial perturbation arises solely from the fluctuations of the curvaton field a 



and neglect the inflaton contribution. Mixed scenario's were considered in ||33|| . After 
the end of inflation the inflaton decays into radiation which dominates the universe. 
We assume the curvaton decays instantaneously into radiation at Hd ec = T, for a 
discussion on the accuracy of this approximation see PH , |55|| . If the curvaton is 
coupled to other scalars, it may also decay non-perturbatively through a parametric 



resonance [36[ [37|, we will not consider this possibility here. 



Using the 5N formalism [39], [40], [41], |42j, the curvature perturbation can be 
expressed in the form 

Ck = N'{t k )Sa k {t k ) + ]-N"{t k ){5a * Sa) k (t k ) + -N"'(t k )(Sa * 6a* 5a) k (t k ) + ... 

2 o 



3 C^osc 3 \ cr nsr \ u, 

/ /// // / \ 

+ r_^ Oosc + 3^050 ) {6(T *s**6(T) k + --- , (2.2) 

where the convolutions are defined by (5a *5a) k (t k ) = (2n)~ 3 J dq5cr q (tfc)<5<7k-q(£fc). 
N(t k ) denotes the number of e-foldings from an initial spatially flat hypersurface t k , 
corresponding to the horizon exit of the mode k, to some final uniform energy density 
surface after the decay of the curvaton. The primes denote derivatives with respect 
to (x(tfc). cr osc sets the scale of the curvaton envelope during the final quadratic 
oscillations before the decay, a(t) = a osc / (mt) 3 ^ 4 , see equation ( |2.7] ) below. r dec 
measures the curvaton contribution to the total energy density at the time of decay, 



3 p. 

'"dec 



A3H 2 



and the results are computed to first order in r^ec throughout this work. 

The information about curvaton self-interactions in (|2.2|) is essentially encoded 
into the derivatives of cr osc . In general, it is not possible to compute cr osc (a(t k )) 
analytically. However, we can obtain some generic information by just looking at the 
evolution equation for the curvaton field in the radiation dominated epoch, 

3 

a + — & + m 2 a + nXa^ 1 = . (2.4) 
2i 



- 3 - 



Switching to the variable x = mt and writing a = a(t k )£(x)x 3//4 , we obtain 

+ ^(1 + ^T 2 ) + f e-^- 3( "- 2)/4 = • (2.5) 
Here we have defined a (relative) self-interaction strength parameter s by 

2\a(t k ) n - 2 



(2.6) 



which is simply the ratio of potential energies stored in the curvaton self-interactions 
and the bare mass part in (|2.1|). Beware that various definitions of s and r dec have 
been used in the literature. 

In the asymptotic limit x — > oo, corresponding to the regime of quadratic oscil- 
lations, ( p.5| ) has a solution £ as = cr osc (s) sin(x + <p{s)), and we obtain the asymptotic 
result 

<r„(t) = a(t k )a osc (s) = crose (mt)3/4 • (2-7) 

From this we learn that <r osc can be expressed in the form 

0osc(o"(4)) = cr(4)<3- osc (s) . (2.8) 



3. Scale-dependent non-Gaussianity 

We analyze the scale-dependence of the non-linearity parameters using the formalism 
developed in [|TJ, |TTJ. The curvaton perturbations at horizon crossing, 5<Tk(t k ), are 
assumed to be Gaussian. We study the effect of relaxing this assumption in Appendix 
We denote the Gaussian part of the curvature perturbation ( |2.2| ) by 

C k ° - N%)Sa k (t h ) = ^ (s) ^2 , (3.1) 

3 a(t k ) 

where the function z(s) is given by 

z(s) = = 1 + — . (3.2) 



The expression (|2.2|) for the curvature perturbation can now be written as 

Ck = C k ° + |/nl(A:)(C G * C G )k + ^9NL(k)(C G * C G * C°)k + • • • , (3.3) 
5 25 

where the non-linearity parameters are given by 

_ 5 N"(t k ) _ 5 / a^ c a OBC \ _ 5 

hM ~ 6 iVW " 4^ I 1 + ^T) = 4^ f{s) ' (3 - 4) 

25 iV"'fa) 25 fCaL^3a>^\ 25 
^ ^ 54 AW = 24^ ^ + "^?- ) = 2l^L » W ' (3 ' 5) 
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to leading order in r<iec- Since we assume ( is generated by a single field, the third non- 
linearity parameter tnl, describing the trispectrum together with #nl, is uniquely 
determined by /nl 

Tnl — 7/nl 



v 5' 

It hence trivially follows that n TNh = 2nf NL , as is always the case for a model where ( 



is generated by a single-source . For a general discussion of the relation between 



the (local) non-linearity parameters see [[43 



Information about the curvaton interactions is encoded into the functions f(s) 
and g(s), which depend only on the self-interaction strength parameter s. For a 
purely quadratic model, s = 0, they read / = 1 and g = since a osc oc cr(tk)- In the 
presence of self-interactions, f(s) and g(s) become non-trivial functions oscillating 
between positive and negative values |22|| . The level of non-Gaussianity may therefore 
strongly deviate from the the naive estimates of |/nl| ~ r dec anc ^ Is'nlI ~ r dcc- 

Equations fl3.4| ) and ( |3.5p are evaluated at the horizon crossing time tk of the 



mode k under consideration, which in general makes /nl and #nl scale-dependent [10 



lTH . The scale dependence can be described by the parameters n/ NL and n 9NL , which 



measure logarithmic derivatives of the non-linearity parameters ([L.2|). Applying the 



results of 1 10 , 111 to the curvaton scenario, we find 



1 N' V' 

i] a f n{n — \){n — 2)s 



f(s) \z(s)(2 + n(n-l)s) 
N" 2 N' V"" 

rj a ( n(n-l)(n-2)s f 3 ^s + n 



g(s) \z(s)(2 + n(n — l)s) \ z(s 

The results are derived to leading order in slow roll, see [|i~0| , jlTf for details. The 
functions z,f,g are defined by ( |3.2| ), (|3.4|) and ( |3.5|) . The slow-roll parameter r\ a 
is defined as usual, r] a = V"/(3H 2 ) = m 2 {2 + n(n — l)s) / '(6H 2 ). In the curvaton 
scenario, the scale dependence of /nl and #nl is entirely generated by curvaton 
interactions. For a purely quadratic model the scale dependence vanishes rif NL = 



n 



9nl 



(note that also r 2 ec g^ = in this case, subleading corrections in r^ec lead 
to #nl = C(/nl)) as the equation of motion for a is fully linear to leading order in 
r dec [0- The first detailed study of the scale dependence of #nl for isocurvature 



models was made in 44 



The results given here are valid for |n/ NL | and |% NL | much less than unity [10, 11]. 
To discuss stronger scale-dependence, the formalism needs to be modified to account 
for the non-Gaussianity of the curvaton perturbations Sa^tj,). A more detailed 
discussion on this issue is presented in Appendix [A| 



- 5 - 



3.1 Regimes of enhanced scale-dependence 

The complicated dynamics of the self-interacting curvaton scenario can lead to a 
considerable enhancement of the scale-dependence. In the regions where \ f{s)\ <C 1 
or \g(s) \ <C 1, the spectral indices 

«/NL « y^y , n 9NL OC (3.8) 

can become much larger than the slow-roll scale r\ a . The amplitudes /nl oc f{s)jr^c 
and #nl oc g(s)/r^ ec , on the other hand, depend not only on the self-interaction 
strength s but also on r dec , measuring the curvaton energy density at the time of its 
decay. As r^cc can be varied independently of s, the non-linearity parameters /nl 
and (?nl can be large even if / or g are suppressed. 

According to [Oj, Planck should be able to probe the scale-dependence of /nl 



to the precision 

50 1 . , 

An /NL ^0.1- =, (3.9) 

where / s k y stands for the fraction of sky observed and the result is derived taking 
/nl — 50,n/ NL = as fiducial values in the analysis. For CMBpol, the error is 
expected to be smaller by a factor of two. Since the error An f NL is inversely propor- 
tional to /nl, it is also interesting to consider the combination /nl^/ nl oc r) a /r<i ec for 
which this dependence drops out. The combination /nl^/ nl is parameterically sup- 
pressed for r\ a <C r dec but can become observable for r) a > r dec . This is in accordance 
with our finding that the regions \f(s)\ <C 1, with f(s) defined in equation ( |3.4| ), 
are characterized by an enhanced scale-dependence. Indeed, the lower bound on r^ec 
following from the observational constraint |/nl| ~ \f(s)\/r dec < 10 2 gets relaxed 



in the regions |/(s)| 1 which makes it easier to have n a > r dcc . For example, if 
\f \ ~ 10 _1 then the observational bound on /nl requires that r^cc ^ 10 -3 , while the 
observed value of the spectral index suggests that r\ a < 1CT 2 , and it could be larger 
if there is an accidental cancellation between r\ a and e. We therefore conclude that 
there is a reasonably large parameter space in which Planck may be able to detect 
both the bispectrum and its scale dependence. 

A similar enhancement of scale-dependence could also take place in other single- 
source models where isocurvature perturbations of an initially subdominant field x 
are converted into curvature perturbations at some later stage, and perturbations of 
the other fields can be neglected. The curvature perturbation can be schematically 
written in a form analogous to the curvaton case J|6|], ( = r(( x + /(s)C^ + •••)> with 
r oc p x / p and s measuring interactions of the x field. This yields /nl oc f(s)/r 
where r, being independent of the time of horizon crossing of a given mode, does not 
contribute to the the scale-dependence. Therefore, ri/ NL = dln|/(s)|/dln/c, and the 
scale-dependence gets enhanced in the regions |/| 1, provided they exist. Similar 
comments apply to the scale-dependence of #nl- 
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3.2 Running of n/ NL and n 9NL 

In the regimes of enhanced scale-dependence, |h./ n J, |n gNL | 3> r) a , the first-order 
derivatives of /nl and <7nl are not necessarily enough to describe the scale-dependence 
but higher-order derivatives may also become relevant. In addition to nj NL and n gNL , 
one then needs to consider the running of these parameters. 

Starting from the expressions (|3.6| ) and (|3.7|) it is straightforward to compute 



the running of n/ NL and n 9NL plj, [44| . The results can be expressed in the form 



= dlnT " + [ 2en ~ 2 + n(n-l)s "*) ^ ' (3 - 10) 

d% NL 2 / 2(n-2)(w-3-3ns) 

d\nk " 9nl+^ £ h ( n _3)(2 + n (n-l)s) 



An 2 (n-2)s f> 
(n-3)(2 + n(n- 1)8)9^^- [6AL) 

where en = —H/H 2 . 

For |n /NL | > 0(rj a ), equation (|3.10|) gives a /NL = -nJ NL + C(e)n /NL , which 
implies that d 2 /NL/dln& 2 = 0(e)df^/d\nk. 0(e) denotes slow roll corrections 
proportional to en or r) a . Similarly all higher order derivatives are slow roll suppressed 
and /nl can be expanded around some reference scale k as 

/nl(*0 = (l + n fNL (k ) ln|- (l + J2 0(O J • (3-12) 

If we consider a range of fc- modes corresponding to a few e- foldings at most, the cor- 
rections O(e)lnk/ko are tiny and can be neglected. This leaves us with the compact 
result 

/nl(*0 = /nl(*o) (l + % NL (^o) hi^J . (3.13) 

Note that this is a non-perturbative expression valid to all orders in n/ NL and not just 
a truncated expansion. A similar result can be derived for g NL whenever |n ffNL | ^> r\ a . 

As a curiosity, we notice that expressions formally similar to ( |3.13| ) appear in 
models where non-Gaussianity is generated by classical superhorizon loops |47j (see 
also [0, f49|). In such scenarios the logarithm \a(k/ko) gets replaced by In (kL) 
where L is an arbitrary infrared cut-off scale |50| . While the role of L is somewhat 



subtle, our expression is manifestly independent of k . 
4. Quartic interactions 

In this Section we discuss curvaton models with (marginally) renormalizable four- 
point interactions 

V =^m 2 a 2 + Xa A . (4.1) 
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In [26] it was shown that for this class of models cr osc can be approximated by, 

i q „— 0.80v / Ao"*/'" i q „— 0.56^ 

ITVn J c . nnr ^, (1 + Ofa,)) , (4.2) 



|r(0.75 + i0.5lVAo-*/m)| |r(0.75 + i 0.36^)1 

where a* = cr(t fe )(l + C(r/ <r )) denotes the curvaton value at the end of inflation. Using 
this result, it is now straightforward to compute the amplitudes /nl and #nl and 
their scale-dependence n/ NL and n 9NL , given by equations (j3.4|)-( ft77| ). The results are 
depicted in Figs. [l] and 0. 




Figure 1: /nl and <?nl plotted against r^ ec , measuring the curvaton energy density at the 
time of decay, and the self-interaction strength parameter s. The contours in the left panel 
run from 10 (black) to 100 (white) with a spacing of 10. In the right panel the contours 
run from —5000 (black) to 1000 (white) with a spacing of 1000; the 0-contours correspond 
to the two horizontal lines. 



|n f ML |/77^ lng NL l/^ 




0.1 0.2 0.5 1 2 5 10 20 0.1 0.2 0.5 1 2 5 10 20 



Figure 2: \nf NL \/r] a and \n gNL \/rj a plotted as a function of the self- interaction strength 
parameter s on logarithmic scales. 



For the quartic case, /nl does not show significant oscillatory features and con- 
sequently rif Nh does not get enhanced for any value of the self-interaction strength 
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parameter s. On the other hand, #nl changes sign twice around s ~ 2, and in this 
region n 9NL is considerably enhanced. The two divergent spikes seen in the plot 
for n gNL correspond to the points where #nl = 0. When moving away from these 
points, #nl starts to grow while n gNL still remains large. In the self-interaction dom- 
inated regime s > 1, both n/ NL and n 9Nh asymptote to constant values. Indeed, in 
this limit equation (O) reduces to a simple power law cr osc oc cr^) 3 / 4 which yields 
f NL = 10/(12r dec ) , n /NL = Ar] a and g NL = -25/(54r| ec ) , n gNL = -20r] a . 

The quartic model nicely demonstrates how the interacting curvaton scenario can 
generate strongly scale-dependent non-Gaussianity. However, the results for n gNL are 
of limited observational interest because the bound |/nl| < 10 2 requires #nl to be 
small |<7nl| ^ 10 3 , see Fig. El This is too small to be detectable with the CMB |51 



5. Non-renormalizable interactions 

For non-renormalizable curvaton potentials, n = 6 and n = 8 in ( |2.1|) , we have used 
numerical methods similar to [^6[ to study the dynamics and compute the scale- 
dependence. In the interaction dominated regime, s ^> 1, it is also possible to obtain 



simple analytical estimates as we briefly discuss at the end of this section, see Sec. 5.1 



The presence of non-renormalizable self-interactions renders both /nl and #nl 



oscillatory functions of the self- interaction strength parameter s pz| . The oscillatory 
behaviour can lead to great enhancement of the scale-dependence, as discussed above. 
This is clearly seen in Figure || which shows \nf NL \/r] a and \n gNL \/r] a as a function of 
s for n = 6. The results for n = 8 are qualitatively similar. 




Figure 3: Plot of \nf NL \/rj a (red-solid line) and \n gNL \/r] a (blue-dashed line) as a function 
of the self-interaction strength parameter s for n = 6. 
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The spikes in the behaviour of |n/ NL | in Figure |3| correspond to points where /nl 
crosses zero. In the vicinity of these points /nl takes non-zero, and for small r^ec 
observable, values while \rif NL \ is one or two orders of magnitude enhanced compared 
to the slow-roll scale r\ a . Similar comments apply to |n SNL | whose behaviour is 
illustrated in the same figure. 

In Figure f| we compare our results with the predicted accuracy of Planck for 
observing the scale-dependence. Figures fJTaTL |4(bj| and |4(c)| show / NL against ^/ NL , 







'r dcc =0.01,tl=0.02 . 
^=0.01.11=0.01 
r^.sO.008, 11=0.005 


i \\ 













^=0.01, r\g= 0.02 
r te =0.01,r,„=0.01 - / 
r Jel .=0.0()8, T| o =0.005 . / 


/ l\ 













(a) n = 4 



(b) n = 6 
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Figure 4: / NL vs. n /NL for n = 4, 6, 8 and g NL vs. n gNL 
r^ec and 77o- . The points shown in the figures range from s 



(d) n = 6 
for ra = 6 with different choices of 



1(T 5 to s = 10 5 . The black 



curves depict the forecasted observational sensitivity of Planck, the region of detectable 
n/ NL lies outside the curves, i.e. further from the origin. 

scanning from s = 10~ 5 to s — 10 5 and keeping r dec and r\ a fixed. Outside the black 
lines which denote |?t-/ nl /nl| = 5, the scale-dependence is detectable by Planck at 
the 1-cr level (and at the 2-a level with CMBPol) fl3|. For n = 6 and n = 8, the 
points corresponding to a given choice of r dec and r\ a do not lie on a single curve as 
in the case n = 4. This is again a manifestation of the oscillatory behaviour of /nl ; 
characteristic for non-renormalizable self-interactions. However, it is noteworthy that 
it is possible to generate observable scale-dependence even for n = 4, provided that 
the ratio "q^jr^ is large enough. 
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For comparison, we have also plotted (?nl against n SNL for n = 6 in Figure |4(d) 



As #nl and rz gNL feel derivatives up to third order, this plot shows considerably more 
structure than the corresponding result for / NL and rij NL (Figure |4(b)| ), which only 
feel derivatives up to second order. There are currently no forecasts on how well n 5NL 
could be measured. 

In Figure |5] we plot n gNL /rif NL against <7nl//nl- These ratios only depend on 
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(c) n = 8 

Figure 5: The observables <7nl//nl an d n 3 N L/ n /NL' which depend on the self-interaction 
strength s only, plotted for s = 10 _5 ...10 5 . The points outside the curves are not accessible 
for any parameter values in the self-interacting curvaton scenario. In the case of a renor- 
malisable self-interaction, n = 4, we also plot points in green calculated using the analytic 
formula. 



the self-interaction strength parameter s. The points depicted in the figure range 
from s = 10~ 5 to s — 10 5 . The predictions asymptote to constant values both 
for s — > and s — > oo and extending the plot region to smaller or larger s-values 
essentially leaves the plots unchanged. Therefore, only the points that lie on the 
curves seen in Figure |5] are accessible in the self-interacting curvaton scenario. The 
region outside the curves can not be accessed for any parameter values. Despite 



its very rich structure and broad range of different observational imprints [22 
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the self-interacting curvaton scenario could therefore be ruled out by a combined 
detection of /nl and #nl and their scale-dependencies. 



5.1 Interaction dominated regime 

In the interaction dominated regime s ^> 1, it is possible to derive analytical results 
even for the non-renormalizable case. In this regime the curvaton oscillations start 
in the non-renormalizable part of the potential and the transition to the quadratic 
potential takes place relatively late after the onset of oscillations. The dynamics can 
therefore be described by the simple scaling law p a oc a ~ 6n ^ n+2 \ unlike for smaller 
values of s p6[j . For s ^> 1, the transition time can be estimated by \a n ~ 2 ~ m 2 
which yields a osc oc cr(tfc)s~ 1 / 8 (l + 0(rj a )). Using this in equations ( pT4| ) - ( ^.7[ ) we 
obtain the results 



1 5(6 -n) Q (rh\ n f _ V* 10(n - 2) / e 2 
r dec 2(10-n) +U {r dc J ' ™ /nl/nL r dec 10 - n + "U 



p2 2(2 -n) , _ fva \ V« 50(6 + n) 



XX 

2 



ta=& siR +0 liJ • "»» 9NL= t« + °feJ • (52) 

The results for n = 4 agree with the discussion in Section ^. For n = 6, the 
amplitudes /nl and g NL vanish to leading order in slow roll, /nl = 0(rj a /r dcc ), 
9nl = fao-Adec)' anc ^ we fi nc ^ h /nl = ^(l)' n 9NL = ^(1)- The first order slow roll 
computation used in our analysis is not enough to derive explicit results for this case 
but the level of non-Gaussianity is clearly unobservably small. For n = 8 we find 
/nl = -5/(2r dec ), #nl = 2/£ L , n /NL = -12^ and n 9NL = 28?7 CT . We have checked 
that these analytical estimates agree with our numerical simulations. 



We note that our results do not agree with those in 29 , in particular compare 



with Eq. (2.56) of [^9| in the large s limit, where the sign of #nl was found to be 
positive for all values of n. We do not attempt to explain the difference, but we note 
again that we have found a good agreement between our analytic and numerical 
results, both here for the interaction dominated regime as well as for the n = 4 case, 
see for example Fig. |5(a) 



6. Conclusions 

We have studied the scale-dependence of the non-linearity parameters, especially 
of /nl and #nl in the curvaton scenario allowing for the possibility of a large self 
interaction. We have found a rich structure in the results, and that a much larger 
scale dependence of the non-linearity parameters is possible than may be expected 
by comparison to the observed spectral index of the power spectrum. This boosts 
the observational prospects for detecting non-Gaussianity and its scale-dependence 
and shows that Planck may achieve a simultaneous detection of /nl and rif Nh . This 
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would put stringent constraints on the curvaton scenario and rule out its simplest 
and most studied version, the curvaton with a quadratic potential. 

Although the richness of the results, as shown in the many plots, makes it hard to 
make firm predictions of the curvaton scenario, one can observe some interesting gen- 
eral trends in the results. In general, increasing the strength of the self interaction s, 
and/or the power of the self-coupling n leads to larger values of the scale-dependence 
as well as an oscillatory structure. Due to their dependence on a third derivative, <?nl 
and n 9Nh have a more complex structure than /nl and n./ NL , which only depend on a 
second derivative. In the limit of no self-interaction, i.e. s = 0, we recover a constant 
and potentially large /nl- However #nl — is far too small to be observable in this 
limit. 

Previous studies of scale-dependence in non-quadratic curvaton scenario's had 
focussed on a region with only small self-interactions. In that regime and for the 



models studied in this article it was found that n/ NL > [2T], while for an axionic 
curvaton potential the opposite sign was found [p^ . In both cases the sign of n/ NL 
was given by the sign of the third derivative of the potential. However we have here 
shown that even for a fixed potential the sign of ny NL may oscillate, depending on 
the initial field value, which affects the self-interaction strength s. For a quartic 
self-interaction there are only oscillations in the sign of n 9NL but not nf NL , for higher 
powers of the self interaction multiple oscillations in both of these parameters occurs. 

The scale dependence of the non-linearity parameters does linearily depend on 
the r\ a slow-roll parameter, just as the power spectrum's spectral index does (provided 
that the e slow-roll parameter is subdominant, one has n s — 1 = 27]^). However the 
numerical coefficient in the case of n/ NL and n 9Nh also depends on the value of the 
self-interaction strength and can become very large in cases where the amplitude of 
the non-linearity parameters become supressed, but they may still remain observable 
provided that the curvaton is sufficiently subdominant at the time of decay. In the 
limit of a large self-interaction strength with an octic self interaction, which does not 
correspond to any suppression of the non-linearity parameters we find n/ NL ~ — 121]^ 
and n 9NL ~ 28r] a , both of which are an order of magnitude larger than the spectral 
index. 

This rich structure of the self-interacting curvaton does not make the model 
unpredictive or unfalsifiable, there are model constraints, for example on how late the 
curvaton decays and on the observed amplitude of the power spectrum which restrict 
the allowed model parameters. The fact that the model can be observationally ruled 
out is clear from the plots relating / NL , gNL, n f NL and n ffNL given in Fig. ||, only a few 
lines in parameter space are allowed. 
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A. On the accuracy of the results 

In computing the scale-dependence, we have neglected the non-Gaussianities of the 
curvaton perturbations Sa^tk) at horizon crossing, following [W, |TTJ. Since we as- 
sume canonical slow-roll dynamics for the curvaton during inflation, the neglected 
parts are in general slow-roll suppressed. (For a discussion of the scale-dependence 
of the three and four-point functions of a test field 5cr k (i fe ), see |53|| .) These non- 
Gaussianities can however play a key role if |n/ N J or |n SNL | become of order unity. 

In this Appendix we address this issue by considering a (unrealistic) toy model 
with a quartic curvaton 

V(a) = Xa A , (A.l) 

and a vanishing classical background field o = 0. The inflationary stage is described 
by a de Sitter solution. While this model cannot lead to a successful curvaton 
scenario, it clearly demonstrates how the non-Gaussianities of 5a^(tk) can become 
important. 

The connected four-point function of the massless curvaton fluctuations, evalu- 
ated at some time t{ after the horizon exit of all the four modes, can be written as 

b m Hi 

(Sa^Sa^Sa^SauM) = (2vr) 3 5(^ k m )^ (7 + £({M) + 

(PikjPik^Pih) + perm.) + 0(X 2 ) . (A.2) 
Here 7 ~ 0.58 is the Euler-Mascheroni constant and £,({k m }) is a dimensionless 



function of the all the four wavenumbers fc m , see [53] for details. P{k) is the spectrum 
of curvaton fluctuations and ki is the mode crossing the horizon at ti. To first order 
in the coupling A, there are no other connected n-point functions (n > 2). 

The four-point function affects the trispectrum of curvature perturbation. Using 
the SN formalism together with ( |A.2j ), we find the non-linearity parameter given 
by 

^N L = |^(l + n° NL ( 7 + e(U + ln^)) • (A.3) 
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Here n° NL denotes the scale-dependence given by (|3 

N' V"" 

9nl jqm 3//2 ' 

(A.4) 

computed assuming the modes 5a^{tk) are Gaussian. (The part proportional to rif NL 
in ( j3.7| ) vanishes as V'"(0) = here.) 

Concentrating, for simplicity, on equilateral configurations k m = k, and setting 
£j = tfe, we obtain 

^ = ^^^1 + ^(7-^ + ^4)^ , (A.5) 
n° 

"»» = l + < L (-,- L f| + l„4) ■ (A - 6) 



The term n° NL (7 - 51/16 + In 4) ~ -1.2 n° NL in (|A.5|) and ( |A.6| ) arises from the 



connected four-point function ( |A.2|) of curvaton fluctuations, that is from the non- 
Gaussianity of <5cr k (t fc ). For |^° NL | <C 1, these corrections can be neglected and we 
recover the results previously used in this work. However, if |n° NL | > 1 the corrections 
clearly have a significant effect on both the amplitude g^L and its scale-dependence. 

It is straightforward to see that the results get modified in a qualitatively similar 
manner for models with realistic potential and a non- vanishing background field a, 

N" n° 
/nl = (1 + 0(4 N J) > % NL = T-^Om . (A.7) 

\ /nl ' 

or Alt// n 

9nl = ^ (1 + ©(njj) , n 9NL = i - . (A.8) 

Therefore, we conclude quite generally that non-Gaussianities of the curvaton per- 
turbations can be safely neglected if In 1 ) I C 1 and |n° I < 1. However, if n° f 

J & 1 /NL 1 1 9NL I ' /NL 

or n° NL become large a more careful analysis is needed, not only to compute the 
scale-dependence but also to find the correct results for /nl and #nl- 
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